Modelling of oedemous limbs and venous ulcers using partial differential equations

BACKGROUND: Oedema, commonly known as tissue swelling, occurs mainly on the leg and the arm. The condition may be associated with a range of causes such as venous diseases, trauma, infection, joint disease and orthopaedic surgery. Oedema is caused by both lymphatic and chronic venous insufficiency, which leads to pooling of blood and fluid in the extremities. This results in swelling, mild redness and scaling of the skin, all of which can culminate in ulceration. METHODS: We present a method to model a wide variety of geometries of limbs affected by oedema and venous ulcers. The shape modelling is based on the PDE method where a set of boundary curves are extracted from 3D scan data and are utilised as boundary conditions to solve a PDE, which provides the geometry of an affected limb. For this work we utilise a mixture of fourth order and sixth order PDEs, the solutions of which enable us to obtain a good representative shape of the limb and associated ulcers in question. RESULTS: A series of examples are discussed demonstrating the capability of the method to produce good representative shapes of limbs by utilising a series of curves extracted from the scan data. In particular we show how the method could be used to model the shape of an arm and a leg with an associated ulcer. CONCLUSION: We show how PDE based shape modelling techniques can be utilised to generate a variety of limb shapes and associated ulcers by means of a series of curves extracted from scan data. We also discuss how the method could be used to manipulate a generic shape of a limb and an associated wound so that the model could be fine-tuned for a particular patient

A survey of partial differential equations in geometric design

YesComputer aided geometric design is an area where the improvement of surface generation techniques is an everlasting demand since faster and more accurate geometric models are required. Traditional methods for generating surfaces were initially mainly based upon interpolation algorithms. Recently, partial differential equations (PDE) were introduced as a valuable tool for geometric modelling since they offer a number of features from which these areas can benefit. This work summarises the uses given to PDE surfaces as a surface generation technique togethe

Multiresolution discrete finite difference masks for rapid solution approximation of the Poisson's equation

YesThe Poisson's equation is an essential entity of applied mathematics for modelling many phenomena of importance. They include the theory of gravitation, electromagnetism, fluid flows and geometric design. In this regard, finding efficient solution methods for the Poisson's equation is a significant problem that requires addressing. In this paper, we show how it is possible to generate approximate solutions of the Poisson's equation subject to various boundary conditions. We make use of the discrete finite difference operator, which, in many ways, is similar to the standard finite difference method for numerically solving partial differential equations. Our approach is based upon the Laplacian averaging operator which, as we show, can be elegantly applied over many folds in a computationally efficient manner to obtain a close approximation to the solution of the equation at hand. We compare our method by way of examples with the solutions arising from the analytic variants as well as the numerical variants of the Poisson's equation subject to a given set of boundary conditions. Thus, we show that our method, though simple to implement yet computationally very efficient, is powerful enough to generate approximate solutions of the Poisson's equation.Supported by the European Union’s Horizon 2020 Programme H2020-MSCA-RISE-2017, under the project PDE-GIR with grant number 778035

Firefly algorithm approach for rational bézier border reconstruction of skin lesions from macroscopic medical images

Image segmentation is a fundamental step for image processing of medical images. One of the most important tasks in this step is border reconstruction, which consists of constructing a border curve separating the organ or tissue of interest from the image background. This problem can be formulated as an optimization problem, where the border curve is computed through data fitting procedures from a collection of data points assumed to lie on the boundary of the object under analysis. However, standard mathematical optimization techniques do not provide satisfactory solutions to this problem. Some recent papers have applied evolutionary computation techniques to tackle this issue. Such works are only focused on the polynomial case, ignoring the more powerful (but also more difficult) case of rational curves. In this paper, we address this problem with rational Bézier curves by applying the firefly algorithm, a popular bio-inspired swarm intelligence technique for optimization. Experimental results on medical images of melanomas show that this method performs well and can be successfully applied to this problem

On the spine of a PDE surface

yesThe spine of an object is an entity that can characterise the object¿s topology and describes the object by a lower dimension. It has an intuitive appeal for supporting geometric modelling operations. The aim of this paper is to show how a spine for a PDE surface can be generated. For the purpose of the work presented here an analytic solution form for the chosen PDE is utilised. It is shown that the spine of the PDE surface is then computed as a by-product of this analytic solution. This paper also discusses how the of a PDE surface can be used to manipulate the shape. The solution technique adopted here caters for periodic surfaces with general boundary conditions allowing the possibility of the spine based shape manipulation for a wide variety of free-form PDE surface shapes

Modelling the Mechanical Behaviour of a Pharmaceutical Tablet Using PDEs.

yesDetailed design of pharmaceutical tablets is essential nowadays in order to produce robust tablets with tailor-made properties. Compressibility and compactibility are the main compaction properties involved in the design and development of solid dosage forms. The data obtained from measured forces and displacements of the punch are normally analysed using the Heckel model to assess the mechanical behaviour of pharmaceutical powders. In this paper, we present a technique for shape modelling of pharmaceutical tablets based on the PDE method. We extended the formulation of the PDE method to a higher dimensional space in order to generate a solid tablet and a cuboid mesh is created to represent the tablet¿s components. We also modelled the displacement components of a compressed PDE- based representation of a tablet by utilising the solution of the axisymmetric boundary value problem for a finite cylinder subject to a uniform axial load. The experimental data and the results obtained from the developed model are shown in Heckel plots and a good agreement is found between both.Available in full text since 5th Feb 2013 following the publisher's embargo period

Cyclic animation using Partial differential Equations

YesThis work presents an efficient and fast method for achieving cyclic animation using Partial Differential Equations (PDEs). The boundary-value nature associ- ated with elliptic PDEs offers a fast analytic solution technique for setting up a framework for this type of animation. The surface of a given character is thus cre- ated from a set of pre-determined curves, which are used as boundary conditions so that a number of PDEs can be solved. Two different approaches to cyclic ani- mation are presented here. The first consists of using attaching the set of curves to a skeletal system hold- ing the animation for cyclic motions linked to a set mathematical expressions, the second one exploits the spine associated with the analytic solution of the PDE as a driving mechanism to achieve cyclic animation, which is also manipulated mathematically. The first of these approaches is implemented within a framework related to cyclic motions inherent to human-like char- acters, whereas the spine-based approach is focused on modelling the undulatory movement observed in fish when swimming. The proposed method is fast and ac- curate. Additionally, the animation can be either used in the PDE-based surface representation of the model or transferred to the original mesh model by means of a point to point map. Thus, the user is offered with the choice of using either of these two animation repre- sentations of the same object, the selection depends on the computing resources such as storage and memory capacity associated with each particular application